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Abstract 

In this article, we consider solutions starting close to some linearly 
stable invariant tori in an analytic Hamiltonian system and we prove 
results of stability for a super-exponentially long interval of time, under 
generic conditions. The proof combines classical Birkhoff normal forms 
and a new method to obtain generic Nekhoroshev's estimates developed 
by the author and L. Niederman in another paper. We will mainly focus 
on the neighbourhood of elliptic fixed points, since with our approach 
the other cases are completely similar. 

1 Introduction and main results 

In this paper, we are interested in the stability properties of some linearly 
stable invariant tori in analytic Hamiltonian systems. Let us begin by the 
case of elliptic fixed points. 

1 . As the problem is local, it is enough to consider a Hamiltonian H defined 
and analytic on an open neighbourhood of in M 2ra , having the origin as 
a fixed point. Up to an irrelevant additive constant and expanding the 
Hamiltonian as a power series at the origin, we can write 

H(z) = H 2 (z) + V(z), 

where z is sufficiently close to in M 2n , Hi is the quadratic part of H at 
and V(z) = 0(\\z\\ 3 ). Recall that the fixed point is said to be elliptic if the 
spectrum of the linearized system is purely imaginary. Then it has the form 
{±ia\, . . . , ±ia n }, for some vector a = (a%, . . . ,a n ) £ W 1 which is called 
the normal (or characteristic) frequency. Due to the symplectic character of 
the equations, such equilibria are the only linearly stable fixed points. Now 
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we assume that the components of a are all distinct so that we can make a 
symplectic linear change of variables that diagonalizes the quadratic part: 

H(z) = f (4 + 4h) + V(z) = "I + V{z), 
i=i 

where / = I(z) is the vector of "formal actions" , that is 

H*) = \(4 + ■ ■ ■ , 4 + 4n) e ^ n 

Assuming the components of a are all of the same sign, it is easy to see that 
H is a Lyapunov function so the fixed point is stable. But in the general 
case, one has to study the influence of the higher order terms V(z), and we 
will explain how it can be done using classical perturbation theory. 

2. Let us first note that, given a solution z(t) of H, if I(t) = I(z(t)) then 

|/(t)|i = £|/i(t)| 
i=i 

is (up to a factor one-half) the square of the Euclidean distance of z(t) to 
the origin, so that Lyapunov stability can be proved if \I(t) — /(0)|i does 
not vary much for all times. 

Now in order to study the dynamics on a small neighbourhood of size 
p around the origin in R 2ra , it is more convenient to change coordinates by 
performing the standard scalings 

z i — > pz, H i — > p~ 2 H, 

to have a Hamiltonian defined on a fixed neighbourhood of zero in M. 2n . 
Then, by analyticity, we extend the resulting Hamiltonian to a holomorphic 
function on some complex neighbourhood of zero in C 2n . So eventually we 
will consider the following setting: we define the Euclidean ball in C 2n 

V s = {ze C 2n | < s} 

of radius s around the origin, and if A s is the space of holomorphic functions 
on V s which are real valued for real arguments, endowed with its usual 
supremum norm | . | s , we consider 

(H(z) = aJ + f(z) 
\H € A s ,\f\s < P- 

Let us emphasize that the small parameter p, which was originally describing 
the size of the neighbourhood of 0, now describes the size of the "perturba- 
tion" / on a neighbourhood of fixed size s. Without loss of generality, we 
may assume s > 3. 
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3. Probably the main tool to investigate stability properties is the con- 
struction of normal forms using averaging methods, and in this case these 
are the so-called Birkhoff normal forms. For an integer m > 1, assuming a 
is non-resonant up to order 2m, that is 

k.a ^ 0, k 6 7U 1 , < \k\i < 2m, 

there exists an analytic symplectic transformation 3> m close to identity such 
that H o <3? m is in Birkhoff normal form up to order 2m, that is 

Ho $ m ( z ) = h m (j) + f m (z), 

where h m is a polynomial of degree at most m in the I variables, and the 
remainder f m is roughly of order p 2 " 1-1 (since before the scaling f m (z) is of 
order ||z|| 2m+1 , see |Bir66] . or |Dou88] for a more recent exposition). The 
polynomials h m are uniquely defined once a is fixed, and are usually called 
the Birkhoff invariants. Therefore the transformed Hamiltonian is the sum 
of an integrable part h m , for which the origin is trivially stable, since I(t) 
is constant for all times, and a smaller perturbation f m . Moreover, if a 
is non-resonant up to any order, we can even define a formal symplectic 
transformation and a formal power series hoc = 5Zfc>i ^ fc > with h m = 
J2k=ih k i such that 

if o$ oc (z) = haa(i). 

In general the series hoc is divergent (this is a result of Siegel) and the 
convergence properties of the transformation are even more subtle (see 
[PM03J). However, Birkhoff normal forms at finite order are still very useful, 
not only because the "perturbation" f m is made smaller, but also because 
the "integrable" part h m , for m > 2, is now non-linear and other classical 
techniques from perturbation theory can be used. 

4. First, in the case n = 2, a complete result of stability follows from KAM 
theory. Indeed, if the frequency a € M 2 is non-resonant up to order 4, the 
Birkhoff normal form reads 

H(z) = a.I + /3I.I + f 2 (z), 

with f3 a symmetric matrix of size n = 2 and f 2 a small perturbation. 
This time we consider the non-linear part /12(f) = 01. 1 + (31. 1 as the inte- 
grable system, and if it is isoenergetically non degenerate, the persistence 
of two-dimensional tori in each energy level close to the fixed point implies 
Lyapunov stability (see [AKN06J, [Arn63], or [Ar n61| for other results). 

However, for n > 3, it is believed that "generic" elliptic fixed points 
are unstable, although this is totally unclear for the moment (see [DLC83J. 
|Dou88j and [KMV04 J). 

Therefore, for n > 3, stability results under general assumptions can only 
concern finite but hopefully long intervals of time, and this is the content 
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of the paper. More precisely, we will prove, under generic assumptions and 
provided p is sufficiently small, that for all initial conditions the variation 
\I(t) — 1(0) |i is of order p for t € T(p), where T{p) is an interval of time 
of order exp (exp(p -1 )) (see Theorem 11.11 for a precise formulation). The 
interpretation in the original coordinates is the following: if a solution starts 
in a sufficiently small neighbourhood of the origin, it stays in some larger 
neighbourhood during an interval of time which is super-exponentially long 
with respect to the inverse of the initial distance to the origin. But first, let 
us describe previously known results on exponential stability, where there 
were basically two strategies. 

5. In a first approach, one assumes a Diophantine condition on a, that is 
there exist 7 > and r > n — 1 such that 

\k.a\ > 7|fc|r, keZ n \{0}, 

but no conditions on the Birkhoff invariants. From the point of view of 
perturbation theory, the linear part is considered as the integrable system. 
In particular, a is non-resonant up to any order, hence we can perform any 
finite number of Birkhoff normalizations, and since we have a control on the 
small divisors, we can precisely estimate the size of the remainder f m (in 
terms of 7 and r). The usual trick is then to optimize the choice of m as 
a function of p in order to obtain an exponentially small remainder with 
respect to p . Therefore the exponential stability is immediately read from 
the normal form, and this requires only an assumption on the linear part (see 
GDF+89] or [DG96]). The ab ove Diophantine condition has full Lebesgue 
measure. However, as we will see later, the threshold of the perturbation and 
the constants of stability are very sensitive to the Diophantine properties of 
a, in particular the small parameter 7. 

6. The second approach is fundamentally different, and it does not rely on 
the arithmetic properties of a. Here, one just assumes that a is non-resonant 
up to order 4, so that the Hamiltonian reduces to 

H(z) = a.I + f3I.I + f 2 (z). 

In this case h,2(I) = a. I + (31. 1 is considered as the integrable system (/3 
being a symmetric matrix of size n). Now we suppose that the non-linear 
part is convex, which is equivalent to (3 being sign definite. Under those 
assumptions, it was predicted and partially proved by Lochak ( |Loc92] and 
[Loc95j ). and completely proved independently by Niederman ( |Nie98| ) and 
Fasso, Guzzo and Benettin ([FGB98J and [GFB98J) that exponential stabil- 
ity holds. Their proofs are based on the implementation of Nekhoroshev's 
estimates in Cartesian coordinates, but they are radically different: the first 
one uses Lochak's method of periodic averagings and simultaneous Diophan- 
tine approximations, while the second one is based on Nekhoroshev's original 
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mechanism. The proof of Niederman was later clarified by Poschel ([P6s99]). 
However, the method of Lochak was restricted to the convex case, and it 
was not clear how to remove this hypothesis to have a result valid in a more 
general context. 

7. In this paper, using the method of [BN09] we are able to replace the 
convexity condition by a generic assumption. Then, combining both Birkhoff 
theory and Nekhoroshev theory as in [MG95J, we will obtain the following 
result. 

Theorem 1.1. Suppose H is as in (HP, with a non-resonant up to any 
order. Then under a generic condition (G) on h^, there exist positive con- 
stants a,a',ci,c<2 and po such that for p < po, every solution z(t) of H with 
|/(0)|i < 1 satisfies 

\I(t) - 1(0) 1 1 < c lP , \t\ < exp (V a ' exp(c 2 a'/r a )) . 

Denoting = Ylk>i ^ anc ^ ^™ = SfcLi ^ ^ us explain our generic 
condition (G) on the formal power series hoc. In fact 

(GQ = |J (G m ) 

m£N* 

consists in countably many conditions, where (G m ) is a condition on h m . 
The first condition (Gi) requires that h\(I) = a. I with a (7, r)-Diophantine 
vector a. The other conditions (G m ), for m > 2, are that each polynomial 
function h m belongs to a special class of functions called SDM^, which was 
introduced in [BN09J (SDM stands for "Simult aneous Diophantine Morse" 
functions, see Appendix [X] for a definition). In this appendix we will show 
that each condition (G m ) is of full Lebesgue measure in the finite dimensional 
space of polynomials of degree m with n variables, assuming r and r' are 
large enough. This is well-known for m = 1, it will be elementary for m = 2 
(see Theorem IA.8|) but for m > 2 it requires the quantitative Morse-Sard 
theory of Yomdin ([Yom83], [YC04], see Theorem IA.3I in the appendix). Let 
us point out that this would have not been possible if we had assumed h m , 
for m > 2, to be steep in the sense of Nekhoroshev, as polynomials are 
generically steep only if their degrees are sufficiently large with respect to 
the number of degrees of freedom (see [LM88j). 

Our condition (G) on the formal series hoc is therefore of "full Lebesgue 
measure at any order". From an abstract point of view, this condition 
defines a prevalent set in the space of formal power series, where prevalence 
is an analog of the notion of full Lebesgue measure in the context of infinite 
dimensional vector spaces. This will be proved in Appendix[Al Theorem lA.6l 
In Theorem II. 1\ we can choose the exponents 

a = (1 + r)-\ a' = 3^(2(n + l)r')- n , 
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and our threshold po depends in particular on 7 and 7'. Moreover our 
constants c\ and C2 also depend on 7 but not on 7', and we shall be a little 
more precise later on. 

As we have already explained, the proof is based on a combination of 
Birkhoff normalizations up to an exponentially small remainder, which are 
well-known (a statement is recalled in Proposition 12. II below), and Nekhoro- 
shev's estimates for a generic integrable Hamiltonian near an elliptic fixed 
point (Theorem [22] below). The latter result is new, and it will follow rather 
easily from the new approach of Nekhoroshev theory in a generic case taken 
in [BN09| . 

8. As a direct consequence of our Nekhoroshev estimates near an elliptic 
fixed point, we can derive an exponential stability result more general than 
those obtained in [FGB98J and |Nie98| . Like in those papers, we only require 
a to be non-resonant up to order 4, and after the scalings 

z 1 — > pz, 77 1 — > p~ 4 H, a 1 — > p 2 a, 

we consider 

H(z)=aJ + f3I.I + f{z) 
HGA s ,\f\s<p. 

However, instead of assuming that /3 is sign definite, our result applies to 
Lebesgue almost all symmetric matrices (3 without any condition on a. Let 
<5 n (R) be the space of symmetric matrices of size n with real entries. 

Theorem 1.2. Suppose 77 is as in IB]) , For Lebesgue almost all (3 E <S n (M), 
there exist positive constants a',b' and po such that, for p < po, every solu- 
tion z{t) of H with |/(0)|i < 1 satisfies 

|7(i) - 7(0) |i < n(n 2 + l)p- b ', \t\ < exp(p- a '). 

The above theorem is a direct consequence of Theorem 12.21 below, pro- 
vided that ^(7) = a. I + pi.I belongs to SDML. But we will prove in 
Appendix [A] that this happens for almost all symmetric matrices /3, inde- 
pendently of a (see Theorem lA.8p . Once again, let us also mention that this 
result is not possible in the steep the quadratic part h 2 (I) = ft 1. 1 is 

steep if and only if f5 is sign definite. In the above statement one can choose 
the exponents 

a' = b' = 3- 1 (2(n + l)rT n , 
and the threshold po depends on 7'. 

9. Let us add that in order to avoid useless expressions, we will only keep 
track of the small parameters p, 7 and 7' and replace any other positive 
constants by a dot (•) when it is convenient. 



6 



Moreover, in this text we shall use various norms for vectors v G M. n 
or v £ C n : | . | will be the supremum norm, | . |i the ^i-norm and | . | the 
Euclidean (or Hermitian) norm. 

10. Let us now describe the plan of the paper. Section [2] is devoted to the 
proof of Theorem 11.11 and Theorem 11.21 In 12.1} we give a statement of the 
Birkhoff normal form up to an exponentially small remainder. In 12.21 we 
will explain how Nekhoroshev's estimates obtained in [B~N0 9j generalize in 
the neighbourhood of elliptic fixed points, and how they imply Theorem ll.2l 
In 12.31 we W1U show how Theorem 11.11 follows from a simple combination of 
Birkhoff 's estimates and Nekhoroshev's estimates, provided our assumption 
on /ioo is satisfied. Then, in section [3l we will state similar results for 
invariant Lagrangian tori and more generally for invariant linearly stable 
isotropic reducible tori. Finally, an appendix is devoted to our genericity 
assumptions. 

2 Proof of Theorem 11.11 and Theorem 11.2 

In the sequel, we recall that we will use the "formal" actions 

i = i(z) = \{4 + ■■■,4 + 4n) e 

but one has to remember that these are nothing but notations for expressions 
in z G M. 2n . Moreover, we will also need to use complex coordinates for the 
normal forms, and, abusing notations, we will also denote them by z € C 2n , 
but of course the solutions we consider are real. 

2.1 Birkhoff 's estimates 

Here we consider a Hamiltonian as in ([2]) , and we assume that the vector a 
is (7, r)-Diophantine. In this context, the following result is classical. 

Proposition 2.1. Under the previous assumptions, if p< 7, then there exist 
an integer m = m(p) and an analytic symplectic transformation 

$ m : £> 3s/4 -> D s 

such that 

Ho<$> m (z) = h m (I) + f m (z) 
is in Birkhoff normal form, with a remainder f m satisfying the estimate 

l/mlss/4 <• />exp (-(7/5 _1 ) a ) , a = (1 + t)-\ 
Moreover, |$ m — Id| 3s /4 <• 7 _1 p and the image of <& m contains the domain 

v s/2 . 
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For a proof, we refer to [GDF+89j and [DG96] . The analo gous result for 



invariant Lagrangian tori can be found in [PW94J or [P6s93j, and in |JV97j 
in the more general case of isotropic and reducible linearly stable invariant 
tori. 

In the above proposition, one has to choose the integer m = m(p) of 
order (7P _1 )^ T+1 - > • So letting p go to zero, the degree of the polynomial 
h m goes to infinity, and this explains why in the proof of Theorem 11,11 we 
require a condition on the whole formal power series h^. 

2.2 Nekhoroshev's estimates and proof of Theorem 11.21 

1. Here we consider the Hamiltonian 

(H(z) = h(I) + f(z) 

[77 € A s , heSDMf, \f\ s <e 

and we have assumed that, on the real part of the domain, the derivatives 
up to order 3 of h are uniformly bounded by some constant M > 1. The 
definition of the set SDML is recalled in Appendix lAl 



Theorem 2.2. Let 77 be as in (E\) , with r' > 2 and 7' < 1. Then there 
exists eo such that if e < Eq, for every solution z{t) with |/(0)| < 1, we have 

\I(t) - /(0)| < (n 2 + l)e b ', \t\ < exp(e- a '), 

with the exponents a' = b' = 3 _1 (2(n + l)r / ) _n . 

Theorem 11,21 is now an immediate consequence of this result and Theo- 
rem IA.8I (see Appendix |Aj) . 

Proof of Theorem From Theorem IA.81 we know that for almost all f3 £ 
S n (R), the Hamiltonian h{I) = a.I + /3I.I belongs to SDM T '(B) with r' > 
n 2 + l. So we can apply Theorem 12 .21 for every solution z(t) with |/(0)| < 1, 
we have 

\I(t) - /(0)| < (n 2 + l)e b \ \t\ < exp(e- a '), 
with the exponents a' = b' = 3 _1 (2(n + l)r')~ n . In particular, this gives 

|7(f) - 7(0)|i < n(n 2 + l)e b ' , \t\ < exp(£- a '), 

for every solution z(t) with |7(0)|i < 1. □ 



2. The statement of Theorem 12.21 is the analogue of the main statement of 
T3N09]. However, the difference is that here we are using Cartesian coordi- 
nates and not action-angle coordinates (i.e. symplectic polar coordinates), 
and we cannot use the latter since they become singular at the origin. So 
we cannot apply directly the main result of [BN09J. This is not a serious 
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issue when applying KAM theory in this context (see [Arn63] or |Pos82j for 
example), but this becomes problematic in Nekhoroshev theory (see |Loc92| 
or [Loc95j for detailed explanations). This result was only conjectured by 
Nekhoroshev in |Nek77| ■ and it took a long time before it could be solved in 
the convex case (|Nie98j,[FGB98j). Here we are able to solve this problem 
in the generic case. The reason is that even though we cannot apply the 
result of [BN09], we can use exactly the same approach, since the method 
of averagings along unperturbed periodic flows is intrinsic, i.e. independent 
of the choice of coordinates, a fact that was first used implicitly in [Nie98j 
and made completely clear in [Pos99| . 

The proof of such estimates usually requires an analytic part, which boils 
down to the construction of suitable normal forms, and a geometric part. 
The geometric part of [BN 09] goes exactly the same way, so in the sequel 
we will restrict ourselves to indicating the very slight modifications in the 
construction of the normal forms. 

3. Consider linearly independent periodic vectors uj\,... ,u n of M. n , with 
periods (T\, . . . , T n ), that is 

Tj = inf{t > | tu>i € Z n }, 1 < j < n. 

Define the domains 

T> r ^ S] (ujj) ={z£ V 8 . I \Vh(I) - Uj\ <■ rj}, l<j<n, 

given two sequences (r±, . . . , r n ) and (s\, . . . , s n ) (recall that V Sj is the com- 
plex ball of radius Sj). We will denote by lj the linear Hamiltonian with 
frequency ujj, that is = ujj.I. 

The supremum norm of a function / defined on T> r ,s (^j) will be simply 
denoted by 

and for a Hamiltonian vector field Xf, we will write 

l^/k^=maxja zi /| r . sr 

To obtain normal forms on these domains we will make the following as- 
sumptions (Aj), for j 6 {1, . . . ,n}, where (A\) is 

{mT\e ■< ri, mT\r\ ■< s\, < r\ <• s\, 
IVi,si(wi) ± 0, si •< s, 

and for j S {2, . . . , n}, (Aj) is 

J mTjE ■< rj, mTjrj ■< Sj, < rj <• Sj, 

\V rjtSj (uj) 0, V r . tS .(u}j) C X> 2rj _ 1 /3,2» J -_ 1 /3(Wi-l). 

With these assumptions, one can prove the following proposition. 



(Ai) 



(A, 
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Proposition 2.3. Consider H = h + / on the domain V ri-Sl (uJi) , with 
\Xf\ ri ,si < £, and let j G {l,...,n}. For any i G {1, if (Ai) is 

satisfied, then there exists an analytic symplectic transformation 

*j : ^2^/3,2^ /a(wj) -> ^..ih) 

ffo^ =h + g j + f j , 
with {gj,k} = /or iG{l,...,j}, and i/ie estimates 

\^gj\2rj/3,2Sj/3 < e i I X fj 1 2^/3,2^/3 <' e 

Moreover, we have $, = $io---ofj 

$i : ^^3,2^/3 (Wi) -4 £V i)Si (w») 

suc/i i/ta* |$* - Id| 2ri / 3i2si /3 < 

The proof is completely analogous to the corresponding one in [BN09], 
Appendix A, to which we refer for more details. In fact, here the proof is 
even simpler since one does not have to use "weighted" norms for vector 
fields. It relies on a finite composition of averagings along the periodic flows 
generated by lj, j G {1, . . . , n}. The case j = 1 is due to Poschel ( [Pos99] ) 
and, for j > 1, the proof goes by induction using our assumption (Aj), 
jG{l,...,n}. 

Once we have this normal form, the rest of the proof in |BN09] goes 
exactly the same way: every solution z(t) of H with |/(0)| < 1 satisfies 

\I{t) - /(0)| < (n 2 + l)e b ', \t\ < exp(e~ a '), 

provided that e < Eq, with £q depending on n, s, M, 7' and r' and with the 
exponents 

a' = b' = 3- 1 (2(n + l)r')" r \ 
2.3 Proof of Theorem fTTTI 

Now we can finally prove Theorem ll.il by using successively Birkhoff's es- 
timates and Nekhoroshev's estimates. 

Proof of Theorem li.il Let H be as in (J2J , first assume that p < p\ with 
pi=-7 so that using our assumption (G?i) we can apply Proposition 12.11 
there exist an integer m = m(p) and an analytic symplectic transformation 

$ m : V 3s/i -> V s 

such that 

H o $ m (z) = h m (I) + f m (z) 
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is in Birkhoff normal form, with a remainder f m satisfying the estimate 

l/m| 3 s/4<-pexp (-(7P -1 ) a ') , a = r + l. 
So let H m = H o $ m , and set 



e = pexp -( 7 p 1 ) a 




By our assumption (G m ), for m > 2, the Hamiltonian // m , which is defined 
on the domain £>3 S /4, satisfies §E\i . Now assume that e <• £o, which gives an- 
other threshold p < p2, with p2 also depending on 7', and our final threshold 
is po = min{/3i, pi\. We can apply Theorem [22} every solution z m (t) of i/ m 
with |J TO (0)|i < 1 satisfies 

\I m (t) -/ m (0)|i <-e 6 ', |t| < exp(e- a '), 

with 

a' = 6' = 3~ 1 (2(n + l)r')~™. 
Recalling the definition of e, this gives 

|/m(t) " An(0)|i <p b ' exp (-6'( 7 p- 1 ) a ) , \t\ < exp (p" ' exp (a'( 7 p- 1 )«)) . 
However, one has 

p b 'exp(-6'(7p- 1 r) <7 _ V 

and as $ m satisfies |$ m — Id| 3s/ / 4 <• 7 _1 p, and its image contains the domain 
I> s /2, a standard argument gives 

- /(0)|! < 1 - 1 p, \t\ < exp exp (a'( 7/3 - 1 ) a )) , 

for any solution z(t) of with |/(0)|i < 1. □ 

3 Further results and comments 

As we have already mentioned, the idea of combining both Birkhoff theory 
and Nekhoroshev theory to obtain super-exponential stability was discovered 
by Morbidelli and Giorgilli ( [MG95] ) in the context of Lagrangian Diophan- 
tine tori. Evidently, we can also state results in this context. 

1. Consider a Hamiltonian system on a manifold which carries an invariant 
Lagrangian Diophantine torus, that is an invariant sub- manifold T which is 
diffeomorphic to the standard torus T n , and whose induced flow is conju- 
gated to a linear flow on T n with a Diophantine frequency. Since the torus 
is Lagrangian, one can locally reduce the situation to a Hamiltonian defined 
011 T*T n = T n x R n , having T n x {0} as the invariant torus. Moreover, by 
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invariance and transitivity of the torus, in the coordinates (9, I) £ T n x R n 
we can write 

H(9,I)=u.I + F(9,I), 

where w is a (7, r)-Diophantine vector and F(9,I) = 0(\I\ 2 ). After some 
scalings one is led to consider 



(C) 



\H(9,I) = u.I + f(6,I) 
[H€A s ,\f\s<P 

where A s is the space of holomorphic functions on the domain 

V s = {(9,1) € (C n /Z n ) x C n I \1{9)\ < s, \I\ < s}, 

with 1(9) the imaginary part of 9. Here one can also define polynomials h ri 
and a formal power series and we can state the following result. 



Theorem 3.1. Suppose H is as in tC]) . Then, under a generic condition 
on hoo, there exist positive constants a, a',ci, c 2 and po such that for p < po, 
every solution (6(t),I(t)) of H with |/(0)| < 1 satisfies 

\I(t) - I(0)| < cip, \t\ < exp \p~ a ' exjp(c 2 a'p^' 



The assumption on /i M and the values of the constants a and a' are the 
same as in Theorem II .1\ as the proof is completely analogous. In fact, it is 
even simpler since we are using action-angle coordinates, and therefore we 
can immediately use Nekhoroshev's estimates obtained in [BN09J without 
any modifications. 

However, it is important to note that one cannot obtain a statement 
similar to Theorem ll.2l simply because in this case a non-resonant condition 
up to a finite order does not allow to build the corresponding Birkhoff normal 
form. 

If we compare this result with [MG95J, our assumption is generic and we 
do not require any convexity. But of course the price to pay is that one has 
to consider the full set of Birkhoff invariants. 

2. As a final result, one can also obtain similar estimates for the general case 
of a linearly stable lower-dimensional torus, under the common assumptions 
of isotropicity and reducibility (which were automatic for a fixed point or a 
Lagrangian torus). In that context, it is enough to consider a Hamiltonian 
defined in T k x R fe x M 21 (by isotropicity), of the form 

H(9, 1, z) = uj.I + \bz.z + F(9, 1, z). 

Here B is a symmetric matrix (constant by reducibility) such that J21B has 
a purely imaginary spectrum (J 2 i being the canonical symplectic structure 
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of M 2i ), and F(9,I,z) = 0(\I\ 2 , ||-z|| 3 )- In those coordinates, the invariant 
torus is simply given by I = 0, z = 0, and this generalizes both the case 
of an elliptic fixed point (where the directions (9, 1) are absent) and of 
a Lagrangian invariant torus (where the directions z are absent). If the 
spectrum {±iai, . . . , ztza/} of J21B is simple, one can assume further that 

H(9, 1, z) = u.I + a.I + F(9, 1, z), 

where I are the "formal actions" associated to the z variables. Therefore, 
after some appropriate scalings we can consider 

H(6,I,z) =u.I + a.I + F(e,I,z) 
HeA s ,\f\ s <P 

where A s is the space of holomorphic functions on the domain 

V s = {(9,1, z) G (C k /Z k ) x C k x C 21 I \1{9)\ < s, \I\ < s, \\z\\ < s}. 

Under a suitable Diophantine condition on the vector (oj,o) G one 
can define polynomials h m and a formal series depending on J = (1,1). 
Birkhoff's exponential estimates in this more difficult situation have been 
obtained in |JV97| . Regarding Nekhoroshev's estimates for a generic inte- 
grable Hamiltonian which depends both on actions and formal actions, they 
can be easily obtained by obvious modifications of our method. Therefore 
we can state the following result. 

Theorem 3.2. Suppose H is as in i fD]) . Then under a generic condition 
on hoo, there exist positive constants a, a', ci, c<i and po such that for p < po, 
every solution (9(t), I(t), z(t)) of H with |J(0)| < 1 satisfies 

\J(t) - J(0)| < cip, \t\ < exp (p~ a ' exp(c 2 aV a )) • 

Once again, the condition on h^ and the values of the exponents are the 
same. 

3. Let us add that one could easily give similar estimates in the dis- 
crete case, that is for exact symplectic diffeomorphisms near an elliptic fixed 
point, an invariant Lagrangian torus or an invariant linearly stable isotropic 
reducible torus. Even if one has the possibility to re-write the proof in 
these settings, the easiest way is to use suspension arguments, as it is done 
qualitatively in |Dou88j or quantitatively in [KP94| (see also |PT97j for a 
different approach) and deduce stability results in the discrete case from the 
corresponding results in the continuous case. 

To conclude, let us mention that important examples of invariant tori 
satisfying our assumptions (linearly stable, reducible, isotropic) are those 
given by KAM theory. However, the latter not only gives individual tori but 
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a whole Cantor family (see [PosOlJ or [AKN06J). In this context, Popov has 
proved exponential stability estimates for the family of Lagrangian KAM 
tori, if the Hamiltonian is analytic or Gevrey ( |PopOO| and |Pop04| ) . His 
proof relies on a KAM theorem with Gevrey smoothness on the parameters 
(in the sense of Whitney) and some kind of simultaneous Birkhoff normal 
form over the Cantor set of tori. We believe that our method should be 
useful in trying to extend those results to obtain super-exponential stability 
under generic conditions. But clearly this is a more difficult problem, and 
the first step is to obtain Nekhoroshev's estimates in Gevrey regularity for 
a generic integrable Hamiltonian, the quasi-convex case having been settled 
in |MS02| . 

A Generic assumptions 

In this appendix, we will show that our assumption (G) is generic, in the 
sense that it defines a prevalent set in the infinite dimensional space of formal 
power series. 

1. Let us first recall the definition of Simultaneous Diophantine Morse func- 
tions (SDM in the following). Let G(n, k) be the set of all vector subspaces 
of W 1 of dimension k. We equip R n with the Euclidean scalar product, 
and given an integer L G N*, we define G L (n,k) as the subset of G(n,k) 
consisting of subspaces whose orthogonal complement can be generated by 
integer vectors with components bounded by L. In the sequel, B will be an 
arbitrary open ball of W l . 

Definition A.l. A smooth function h : B — > 1R is said to be SDM if there 
exist 7' > and r' > such that for any L G N*, any k G {1, . . . ,n} 
and any A G G L (n,k), there exists (ei,...,e&) (resp. (/1, . . . , f n -k))> an 
orthonormal basis of A (resp. of A- 1 ), such that the function defined on 
B by 

h A (a, (3) = h (aiei H h a k e k + H h Pn-kfn-k) 

satisfies the following: for any G B, 

\\d a h&(a, < j'L~ T ' =^ \\d aa h A (a,P).r]\\ > i L~ T ' ||r/||, 
for any n G R h \ {0}. 

This definition is inspired by the steepness condition of Nekhoroshev and 
the quantitative Morse-Sard theory of Yomdin (see [BN09j for more expla- 
nations). It depends on a choice of coordinates adapted to the orthogonal 
decomposition A® A- 1 , so for A G G L (n,k) and (a, /3) G B, d a h A (a,(3) is 
a vector in M. k and d aa h\{a, (3) is a symmetric matrix of size k with real 
entries. 
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Remark A. 2. Note also that the definition can be stated as the following 
alternative: for any (a, (3) £ B, either we have \\d a h\(a, (3)\\ > 7-L~ r or 
\\daahx{oi, > 7-^~ T |MI for any rj € M. k \{0}. Hence for a given function 

it is sufficient to verify that \\d aa h^{a, /3).?y|| > 7L _T ||7y|| for any rj E IR fc \{0} ; 
and we will use this fact later (in Theorem \A.8\) . 

2. The set of SDM functions on B with respect to 7' > and r' > will 
be denoted by SDM^, (B), and we will also use the notation 

SDM T '{B) = [j SDM^(B). 

7'>0 

The following theorem was proved in [BN09J, and it relies on non trivial 
results from quantitative Morse-Sard theory ( [Yom83] . [YC 04J ) . 

Proposition A. 3 f |BN09| l. Let r > 2(n 2 + 1), and h G C 2n+2 (B). Then 
for Lebesgue almost all £ 6 IR n , the function h^, defined by h^(I) = h{I) — ^.I 
for I £ B, belongs to SDM T '(B). 

Now let us recall the definition of a prevalent set ([HSY92J, see also 
[UY05Q . 

Definition A. 4. Let E be a completely metrizable topological vector space. 
A Borel subset S C E is said to be shy if there exists a Borel measure \x on 
E, with < /i(C) < 00 for some compact set C C E, such that fi(x + S) = 
for all x £ E. 

An arbitrary set is called shy if it is contained in a shy Borel subset, and 
finally the complement of a shy set is called prevalent. 

For a finite dimensional vector space E, by an easy application of Fubini 
theorem, prevalence is equivalent to full Lebesgue measure. The following 
"genericity" properties can be checked ( [OY05] ) : a prevalent set is dense, a 
set containing a prevalent set is also prevalent, and prevalent sets are stable 
under translation and countable intersection. Furthermore, we have an easy 
but useful criterion for a set to be prevalent. 

Proposition A. 5 ([HSY92]). Let A be a Borel subset of E. Suppose there 
exists a finite dimensional subspace F of E such that, denoting Xf the 
Lebesgue measure supported on F, the set x + A has full Xf -measure for 
all x € E. Then A is prevalent. 

It is an obvious consequence of Proposition [AT3l and Proposition IA.5I that 
SDM T '(B) is prevalent in C 2n+2 (B) for r' > 2(n 2 + 1). 

3. Now let Voo = M[[-Xi, . . . ,X n ]\ be the space of all formal power series 
in n variables with real coefficients. It is naturally a Frechet space, as the 
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projective limit of the finite dimensional spaces V m consisting of polynomials 
in n variables of degree less than or equal to m. We define the subset 

S T ^ = {h^ G Voo | h m G SDM T '(B), Vm > 2}, 

where h m = YlT=i h k if ^oo = Sfe>i h k , and we identify the polynomial h m 
with the associated function defined on B. Let us also define 

= {h^ G Voo | h!(X) = a.X, a G V T }, 

where T> T is the set of Diophantine vectors of M" with exponent r, and finally 

££> T = v lo n5^. 

The set Q T ^ is the set of formal power series for which condition (G) holds. 

Theorem A. 6. For r > n — 1 and t' > 2(n 2 + 1), the set Goo is prevalent 
in Voo- 

Proof. As the intersection of two prevalent sets is prevalent, it is enough to 
prove that both sets V^, for r > n — 1, and <S^o, for r' > 2(n 2 + 1), are 
prevalent. 

For the set , this is an easy consequence of the fact that D T is of full 
Lebesgue measure in M n , for r > n— 1, and Proposition IA.5I with F = Vi, 
the space of linear forms. For the set S^, first note that we can write 



^oo I I ^oo,m> 



m>2 

where, for an integer m > 2, 

«S£,m = {^oo G Poo I h m G SL>M T '(B)}. 

As a countable intersection of prevalent sets is prevalent, it is enough to 
prove that for each m > 2, the set S^ om is prevalent in Poo- But once 
again this is just a consequence of Proposition IA.3I and Proposition IA.5I 
with F = V\ the space of linear forms. □ 

For m > 2, the set of polynomials h m for which condition (G m ) is satis- 
fied is given by 

= {h m £V m \h m £ SDM T '(B)}, 
and the proof of the above theorem immediately gives the following result. 
Theorem A. 7. For r' > 2(n 2 + 1), the set is of full Lebesgue measure 

in l^rn . 
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4. Now in the special case m = 2, we can state a refined result which is 
due to Niederman ( |Nie07j ) . 

Theorem A. 8. For Lebesgue almost all (3 £ «S n (K), the function 

h(I) = a.I + pi.I 

belongs to SDM T ' (B) provided r' > n 2 + 1. 

In the above theorem, there is no condition on a, and contrary to Propo- 
sition IA.31 the proof does not rely on Morse-Sard theory. Let us denote by 
A the one-dimensional Lebesgue measure and by 1^ the identity matrix of 
size k. We shall use the following elementary lemma. 

Lemma A. 9. Let k £ {1, . . . ,n}, (3^ £ Sfc(R) and k > 0. Then there exists 
a subset C K CR such that 

\{C K ) < 2kn, 

and for any £ ^ C K , the matrix (3^^ = Pk ~ £,Ik satisfies 

Wk£-V\\ > 

for any rj £ R k \ {0}. 

Of course, the set C K depends on the matrix 

Proof. Let {Ai,... , A&} be the eigenvalues of then in an orthonormal 
basis of eigenvectors for (3k, the matrix (3k£ is also diagonal, with eigenvalues 
{Ai — Afc — £}. Then one has HA-^.^H > k\\i]\\ for any r\ £ M k \ {0} 

provided that for all i £ {1, . . . , k}, |Aj — £| > k, that is if £ does not belong 
to 



C K = |J[Aj - K,\i + «]. 



The measure estimate A(C K ) < is trivial. □ 

With this lemma, the proof is now similar to that of Proposition IA.3I 

Proof of Theorem [A~R Let h(I) = a. I + (31. 1, and given A £ G L (n,k), 
we denote by /3a £ Sfc(lR) the matrix which represents the quadratic form 
(31. 1 restricted to the subspace A. Since the second derivative of h along 
any subspace is constant, then coming back to definition IA.1I and using 
remark \EM h £ SDMjf if 

\\P A . V \\ >>y'L- r '\\r,\\, (1) 

for any A £ G L (n, k) and any rj £ R k \ {0}. Let A T ', be the subset of S n (R) 
whose elements contradict condition (P), that is 

A^,={/3ES n (R) | \\(3 A .r,\\ < j'L^' \\r,\\, AeG L (n,k), r, £ R k \ {0}}, 
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and 

a* = n <. 

7'>0 

What we need to show is that A T has zero Lebesg ue measure in S n 
provided r' > n 2 + 1. 

Apply lemma DOJl to /3a G Sfc(R), with k = r ]'L~ T , to have a subset 
C 7 ', r ',L,A C R such that 

A(Cy ir / )L) A) < 2fc 7 'i" T ', (2) 
and for any £ ^ Cy^i^A) the matrix /3a,£ = /3a — satisfies 

H^.7711 >7^- T '||r?ll 
for any i)£l l "\ {0}. If we define 

= U U U ^7',t',L,A' 

LeN* fce{l,...,n} AeG L (n,fc) 

then 

C 7 ,, r , = {£ e R | = /3 - £In e 

and so 

C T / = P| Cy jT * = {£ G R | ^ G v4 T '}. 

7'>0 

It remains to prove that the Lebesgue measure of C T / is zero, since by Fubini 
theorem, this will imply that the Lebesgue measure of A T is zero. By our 
estimate ([2]), we have 

n 

a(c 7V o < J2\ GL ( n >w 2k 7' L ~ T ' 

LeN* fc=i 

n 

LeN* fc=i 

U=i / \LeN* / 



n(n + 1) ( £ L« 2 -^ y 

\LeN* / 



and, since r' > n 2 + 1, the above series is convergent. Hence 

A(C r /) = inf A(Cy r ') = 0. 

7'>0 



□ 
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